Double giant resonances in deformed nuclei 
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to 1.33 of that for the single resonance. 
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One of the most exciting progresses in the field of gi- 
ant resonances in atomic nuclei for the last few years was 
the experimental observation of two-phonon giant reso- 
nances M. Nowadays, we may speak about some sys- 
tematics of their properties (the energy position, width 
and excitation probability) in spherical nuclei although 
it is still sparse and some open questions in this field 
stimulates theoretical studies (see, e.g. Refs. |p|4l"o||). In- 
vestigation on the properties of two-phonon giant reso- 
nances together with similar studies on low-lying two- 
phonon states 1 1 1 1 should give an answer on how far the 
harmonic picture of boson-type excitations holds in the 
finite fermion systems like atomic nuclei. 

The possibility to observe two-phonon giant resonances 
in deformed nuclei with the present state of art experi- 
mental techniques is still questionable. This is mainly 
due to the fact that one has to expect a larger width of 
these resonances as compared to spherical nuclei. Also, 
the situation with the low-lying two-phonon states in de- 
formed nuclei is much less clear than in spherical ones. 

The first experiment with the aim to observe the dou- 
ble giant dipole resonance (DGDR) in 238 U in relativistic 
heavy ion collisions (RHIC) was performed recently at 
the GSI/SIS facility by the LAND collaboration @. It 
will take some time to analyze the experimental data and 
to present the first experimental evidence of the DGDR 
in deformed nuclei, if any. The first microscopic study of 
properties of the DGDR in deformed nuclei is the subject 
of the present paper. The main attention will be paid to 
the width of the DGDR and its shape. 

In a phenomenological approach the GDR is considered 
as a collective vibration of protons against neutrons. In 
spherical nuclei this state is degenerate in energy for dif- 
ferent values of the spin J = 1~ projection M = 0, ±1. 
The same is true for the 2 + component of the DGDR 
with projection M = 0,±1,±2. In deformed nuclei with 
an axial symmetry like 238 U, the GDR is spit into two 
components I n (K) = 1~(0) and I n (K) = 1"(±1) cor- 
responding to vibrations against two different axes. In 
this approach one expects a three-bump structure for the 
DGDR with the value K = 0, K = ±1 and K = 0,±2, 
respectively. Actually, the GDR possesses a width and 
the main mechanism responsible for it in deformed nu- 
clei is the Landau damping. Thus, the conclusion on 
how three bumps overlap and what is the real shape of 
the DGDR in these nuclei, i.e., either a three-bump or a 
flat broad structure, can be drawn out only from some 
consistent microscopic studies. 

In the present paper we use the Quasiparticle Phonon 
Model (QPM) |13) to investigate the properties of the 
GDR and the DGDR in 238 U. The QPM, although with 
somewhat different technical details which reflect the dif- 
ference between spherical and deformed nuclei, was used 
to investigate the same resonances in spherical nuclei in 
Refs. ^.^,||,|o|). The model Hamiltonian includes an av- 
erage field for protons and neutrons, monopolc pairing 
and residual interaction in a separable form. We use 
in our calculations for 238 U the parameters of Woods- 



Saxon potential for the average field and monopole pair- 
ing from the previous studies They were adjusted 
to reproduce the properties of the ground state and low- 
lying excited states. The average field has a static defor- 
mation with the deformation parameters (3% = 0.22 and 
/?4 = 0.08. To construct the phonon basis for the K = 
and K — ±1 components of the GDR we use the dipolc- 
dipole residual interaction (for more details, see e.g. 
Ref. [^3|). The strength parameters of this interaction 
are taken from Ref. |l5| where they have been fitted to 
obtain the centroid of the B(E1, 0+ s -> l~ (K = 0,±1)) 
strength distribution at the value known from experiment 
p6| and to exclude the centre of mass motion. In this 
approach, the information on the phonon basis (i.e. the 
excitation energies of phonons and their internal fermion 
structure) is obtained by solving the RPA equations. For 
electromagnetic El-transitions we use the free values of 
the effective charges, %ff ~ eN(—Z)/A. 

The results of our calculation of the B(E1) strength 
distribution over |lx_a(i)) an d |l^ = ± 1 («')) GDR states 
are presented in Fig. |l|, together with the experimental 
data. The index i in the wave function stands for the 
different RPA states. All one-phonon states with the en- 
ergy lower than 20 MeV and with the B(E1) value larger 
than 10 -4 e 2 /m 2 are accounted for. Their total number 
equals to 447 and 835 for the K = and K = ±1 com- 
ponents, respectively. Only the strongest of them with 
B(E1)> 0.2 e 2 fm 2 are shown in the figure by vertical 
lines. Our phonon basis exhausts 32.6% and 76.3% of 
the energy weighted sum rules, 14.8 • NZ/A e 2 fm 2 MeV, 
by the K = and K = ±1 components, respectively. For 
a better visual appearance we also present in the same 
figure the strength functions averaged with a smearing 
parameter, which we take as 1 MeV. The long (short) 
dashed-curve represent the K = (K = ±1) components 
of the GDR. The solid curve is their sum. The calcula- 
tion reproduces well the two-bump structure of the GDR 
and the larger width of its K = ±1 component. The last 
is consistent with the experiment jl6| which is best fitted 
by two Lorentzians with widths equal to Li = 2.99 MeV 
and T2 = 5.10 MeV, respectively. The amplitudes of both 
maxima in the calculation are somewhat overestimated 
as compared to the experimental data. This happens 
because the coupling of one-phonon states to complex 
configurations is not taken into account which can be 
more relevant for the K — ±1 peak at higher energies. 
But in general the coupling matrix elements are much 
weaker in deformed nuclei as compared to spherical ones 
and the Landau damping describes the GDR width on a 
reasonable level. 

The wave function of the + and 2 + states belong- 
ing to the DGDR are constructed by the folding of two 
1~ phonons from the previous calculation. When a two- 
phonon state is constructed as the product of two iden- 
tical phonons its wave function gets an additional factor 
\j\f2. The 1 + component of the DGDR is not consid- 
ered here since its excitation is quenched in RHIC for 
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the same reasons as in spherical nuclei M. The anhar- 
monicity effects which arise from interactions between 
different two-phonon states are also not included in the 
present study. It was shown that these effects have an 
yj-4/3 dependence on the mass number A M and that 
they are small for the DGDR @,|J| even for 136 Xe and 

208p b 

The folding procedure yields three groups of the 
DGDR states: 



a) [ljf=o(*i)® 1 k=o(*2)] + =0 , 2 + =0 ) 

b) |[lx=o(i)®l*=±i(0]a+ =±1 ) 

C) |[lK= ± l(il)®lK= ± l(4)] + =0 ,2+ =Di±2 )- (0-1) 



and 



The total number of non-degenerate two-phonon states 
equals to about 1.5 • 10 6 . The energy centroid of the first 
group is the lowest and of the last group is the highest 
among them. So, we also obtain the three-bump struc- 
ture of the DGDR. But the total strength of each bump is 
fragmented over a wide energy region and they strongly 
overlap. 

Making use of the nuclear structure elements discussed 
above, we have calculated the excitation of the DGDR 
in 238 U projectiles (0.5 GeV-A) incident on 120 Sn and 
208 Pb targets, following the conditions of the experiment 
in Ref. |L2 j. These calculations have been performed in 
the second order pert urb ation theory [Q, in which the 
DGDR states of Eq. (3.1) are excited within a two-step 
process: g.s.^GDR^DGDR. As intermediate states, the 
full set of one-phonon 1 1^ =0 (*)) an d |lif=±i(*')) states 
was used. We have also calculated the GDR excitation 
to first order for the same systems. The minimal value 
of the impact parameter, which is very essential for the 
absolute values of excitation cross section has been taken 



according to 6„ 



i.28- (a] /3 + a; /0 ) 



1/3, 



The results of our calculations are summarized in Fig. |2j 
and Table ffl. In Fig. || we present the cross sections of 
the GDR (part a) and the DGDR (part b) excitation in 
the 238 U (0.5 GeV-A) + 208 Pb reaction. We plot only 
the smeared strength functions of the energy distribu- 
tions because the number of two-phonon states involved 
is numerous. The results for 238 U (0.5 GeV-A) + 120 Sn 
reaction look very similar and differ only by the abso- 
lute value of cross sections. In Table | the properties of 
the GDR and the DGDR, and their different K compo- 
nents are given. The energy centroid E c and the second 
moment, m 2 = \ZJ2k a k ' Wk ~ E c) 2 1 J2k a k, of the dis- 
tributions are averaged values for the two reactions under 
consideration. 

The two-bump structure can still be seen in the curve 
representing the cross section of the GDR excitation in 
238 U in RHIC as a function of the excitation energy. But 
its shape differs appreciably from the B(E1) strength dis- 
tribution (see Fig. Ha in comparison with Fig. g). The 
reason for that is the role of the virtual photon spectra. 



First, for the given value of the excitation energy and 
impact parameter it is larger for the K — ±1 compo- 
nent than that for the K = one (see also the first two 
lines in Table |). Second, for both components it has 
a decreasing tendency with an increase of the excitation 
energy |l7). As a result, the energy centroid of the GDR 
excitation in RHIC shifts by the value 0.7 MeV to lower 
energies as compared to the same value for the B(E1) 
strength distribution. The second moment mi increases 
by 0.2 MeV. 

The curves representing the cross sections of the ex- 
citation of the K = ±1 and K = ±2 components of 
the DGDR in 238 U in RHIC have typically a one-bump 
structure (see the curves with squares and triangles in 
Fig. ||b, respectively). It is because they are made of two- 
phonon 2 + states of one type: the states of Eq. ([Tljb) and 
Eq. (|0.l| c), respectively. Their centroids should be sepa- 
rated by an energy approximately equal to the difference 
between the energy centroids of the K = and K = ±1 
components of the GDR. They correspond to the second 
and the third bumps in a phenomenological treatment of 
the DGDR. The K = components of the DGDR include 
two group of states: the states represented by Eq. (0.1a) 
and those of Eq. (0.1c). Its strength distribution has two- 
bumps (see the curve with circles for the 2 + (K = 0) and 
the dashed curve for the + (K = 0) components of the 
DGDR, respectively). The excitation of the states given 
by Eq. (|0.l| a) in RHIC is enhanced due to their lower 
energies, while the enhancement of the excitation of the 



states given by Eq. (3.1c) is related to the strongest re- 
sponse of the K = ±1 components to the external El 
Coulomb field in both stages of the two-step process. 

Summing together all components of the DGDR yields 
a broad one-bump distribution for the cross section for 
the excitation of the DGDR in 238 U, as a function of exci- 
tation energy. It is presented by the solid curve in Fig. ||b. 
Another interesting result of our calculations is related to 
the position of the DGDR energy centroid and to the sec- 
ond moment of the DGDR cross section. The centroid of 
the DGDR in RHIC is shifted to the higher energies by 
about 0.8 MeV from the expected value of two times the 
energy of the GDR centroid. The origin for this shift is 
in the energy dependence of the virtual photon spectra 
and it has nothing to do with anharmonicities of the two- 
phonon DGDR states. In fact, the energy centroid of the 
B(E1, g.s. -> lr)x B(E1, ir — >DGDR/) strength func- 
tion appears exactly at twice the energy of the centroid 
of the B(E1, g.s. -^GDR) strength distribution because 
the coupling between different two-phonon DGDR states 
are not accounted for in the present calculation. The 
same shift of the DGDR from twice the energy position 
of the GDR in RHIC also takes place in spherical nu- 
clei. But the value of the shift is smaller there because in 
spherical nuclei the GDR and the DGDR strength is less 
fragmented over their doorway states due to the Landau 
damping. For example, this shift equals to 0.25 MeV in 
208 Pb for the similar reaction. This effect is also seen 
when the DGDR position against the GDR is reported 
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from experimental studies |18|]. But the larger value of 
the shift under consideration in deformed nuclei should 
somehow simplify the separation of the DGDR from the 
total cross section in RHIC. 

Another effect which also works in favor of the extrac- 
tion of the DGDR from RHIC excitation studies with 
deformed nuclei is its smaller width than times the 
width of the GDR, as observed with spherical nuclei. Our 
calculation yields the value 1.33 for the ratio rDGDR/roDR 
in this reaction. The origin for this effect is in the dif- 
ferent contributions of the GDR K = and K = ±1 
components to the total cross section, due to the reac- 
tion mechanism. It should be remembered that only the 
Landau damping is accounted for the width of both the 
GDR and the DGDR. But we think that the effect of 
narrowing of the DGDR width still holds if the coupling 
to complex configurations is included in the calculation. 

To conclude, we present the first theoretical studies 
based on microscopic calculation of the properties of the 
two-phonon giant dipole resonance in deformed nuclei in 
relativistic heavy ion collisions. We predict that the ex- 
citation function has a one-bump shape and that there 
are at least two effects which work in favor of its exper- 
imental observation, namely, the energy shift to higher 
energies, and the narrowing of its width. 
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TABLE I. The properties of the different components of 
the GDR and the DGDR in 238 U. The energy centroid E c , 
the second moment of the strength distribution m 2 in RHIC, 
and the cross sections a for the excitation of the projectile 
are presented for: a) 238 U (0.5 GeV-A) + 120 Sn, and b) 238 U 
(0.5 GeV-A) + 208 Pb. 







E c 
[MeV] 


m 2 
[MeV] 


a [mb] 
a) b) 


GDR(_R- = 0) 




11.0 


2.1 


431.2 


1035.4 


GDR(K = ±1] 


i 


12.3 


2.6 


1560.2 


3579.1 


GDR(total) 




12.0 


2.6 


1991.4 


4614.5 


DGDR +(K = 


0) 


25.0 


3.4 


18.3 


88.9 


DGDR 2+ (if = 


0) 


24.4 


3.5 


11.8 


58.7 


DGDR 2 +(if = 


±1) 


23.9 


3.2 


22.7 


115.4 


DGDR 2+ (K = 


±2) 


25.3 


3.4 


49.7 


238.3 


DGDR(total) 




24.8 


3.4 


102.5 


501.3 
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FIG. 1. The B(E1) strength distribution over K = 
(short-dashed curve) and K = ±1 (long-dashed curve) 1~ 
states in 238 U. The solid curve is their sum. The strongest 
one-phonon 1~ states are shown by vertical lines, the ones 
with K — are marked by a triangle on top. Experimental 
data are from Ref. flirt. 



FIG. 2. The strength functions for the excitation: a) of the 
GDR, and b) of the DGDR in 238 U in the 238 U (0.5 GeV-A) 
_l_ 208P)-, reaction. In a), the short-dashed curve corre- 
sponds to the GDR (K — 0) and the long-dashed curve to 
the GDR (K — ±1). In b) the dashed curve corresponds 
to the DGDR + (K = 0), the curve with circles to the 
DGDR 2 +(K = 0), the curve with squares to the DGDR 2 + 
(K — ±1), and the curve with triangles to the DGDR 2 + 
(K — ±2) . The solid curve is the sum of all components. 
The strength functions are calculated with the smearing pa- 
rameter equal to 1 MeV. 



